A FOUR-DIMENSIONAL HOOKE'S LAW CAN 
ENCOMPASS LINEAR ELASTICITY AND INERTIA 
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O^ . Abstract. The question is examined, whether the formaUy straight- 






forward extension of Hooke's time-honoured stress-strain relation to the 
four dimensions of special and of general relativity can make physical 
sense. The four-dimensional Hooke's law is found able to account for 
C ' the inertia of matter; in the flat space, slow motion approximation the 

^ , field equations for the "displacement" four-vector field ^* can encom- 

pass both linear elasticity and inertia. In this limit one just recovers the 
CO ' equations of motion of the classical theory of elasticity. 

(N 

>, 

^ , 1. Introduction 

Os . 

^^ ' After having fostered the birth of the four-dimensional approach to special 

relativity ^ , macroscopic electromagnetism has found its natural expression 

Q\ \ in the four dimensional language of general relativity, both in vacuo [Q and 

0^ ■ in matter S, Q. Its field quantities and its field equations have achieved 

O , a canonical form, that can be summarized as follows 0: the unconnected 

^^1 space-time manifold suffices for writing Maxwell's equation in the naturally 

5h . invariant form: 



X 



(1.1) u% = s^ 

(1-2) F"[ik,m] = 0: 

where H*'^ is a skew, contravariant tensor density that represents the electric 
displacement and the magnetic field, while -Fj^ is a covariant skew tensor 
that accounts for the electric field and for the magnetic induction. Even if 
the four-current density s* is given and the co-ordinate system is fixed, these 
equations are not sufficient for determining both H and -Fjfc. They fulfil the 
two identities H^\ i ~ ^^ which ensures the conservation of the electric four- 
current, and e "Kjfc „ji „ = 0, where e " is the totally antisymmetric 
tensor density of Ricci and Levi Civita. By simply counting the components 
of the fields, the equations and the identities one gathers that ( |l.l| ) and 



(|l.2D need to be complemented by the so-called constitutive equation, i. e. 
by some tensor equation that uniquely defines for instance H in terms of 
Fik and of whatever fields may be needed in describing the features of the 
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electromagnetic medium. For a linear medium the constitutive equations 
can be written as |Q: 

(1.3) n'" = ix^'^-^F^^; 

the properties of the medium are specified by the four-index tensor den- 
sity X*'^™". To complete this schematic picture of macroscopic electromag- 
netism, one may add that, once the constitutive equation is given, Maxwell's 
equations can be solved in terms of the four- vector potential ipi. If Fj^ is 
defined as: 

(1-4) Fik = fk,i - ^i,k, 



equation ( |1.2| ) is satisfied. Then the constitutive equation (^) takes the 
form: 

(1.5) ff ^ = ^X"^"'^^n,m - ^m,n) 

and the four components of (/?j are determined, up to a gauge transformation. 



by the four equations (1.1). One may well wonder why the introduction to 
a paper announced to deal with elasticity and inertia begins with these elec- 
tromagnetic reminiscences. But equation (|1.5D , this pervasive, general rel- 
ativistic, four-dimensional constitutive equation is the antisymmetric coun- 
terpart of an equally pervasive, although pre-relativistic, three-dimensional 
and symmetric constitutive equation: 

(1-6) Q''' = lc'^'''{^p,. + &,p)- 

This is the law of linear elasticity, disclosed |^] by Hooke in 1678 with the 
words "iti tensio sic vii\ It is given here by using the language of three- 
dimensional tensors and of tensor analysis, that were just invented to cope 
with the far-reaching developments stemmed from Hooke's discovery. In 
this way all the analogies between (|1.5| ) and ( |1.6| ) are made apparent: the 
four-potential (pi{x^) is the four-dimensional counterpart of the displacement 
vector ^p{x ), the deformation tensor 



1 



(1-7) V=o(W + ^- 



2 



-,pJ 



is the symmetric, three-dimensional counterpart of Fik, while the symmet- 
ric, three-dimensional stress tensor 0'*'^ of (1.6) is replaced in ( |1.5| ) with the 



skew, four-dimensional tensor density H*'^. One might indulge in disclosing 
further analogies that betray the mechanistic origin of electromagnetism, but 
a question comes to the mind: have equation ( |1.6| ) and the classical theory 
of elasticity found a generally accepted reformulation in the framework of 
special and of general relativity, as it has occurred with electromagnetism? 
Even a cursory inspection of the literature shows that the answer is negative. 
This problem has attracted considerable attention in the years soon after 
1905, as it is testified e. g. by the seminal paper of Herglotz on the mechan- 
ics of deformable bodies Q. After 1915 Nordstrom resumed |^] Herglotz' 



approach and translated it in a general relativistic form, but, to our knowl- 
edge, the interest in the issue sinked, and only resurfaced at the time when 
Weber undertook Q his studies on the detection of gravitational waves. In 
classical elasticity, the deformation of an elastic body is measured relative 
to a natural unstrained state. Synge found it difficult to carry over this 
idea into a pseudo-Riemannian space-time and decided Q that the sentence 
"rate of change of stress linear function of rate of strain" should replace, in 
general relativity, Hooke's Latin dictum. Rayner resurrected the concept of 
reference state under the form of a reference metric [10| and postulated the 
following relativistic variant of Hooke's law: 



(1.8) 




s.,^ 


^Cifi-fki-l 


where the auxiliary 


metric 






(1.9) 




y^ 


= g'^ + u'vf' 



defined in terms of the true metric git and of the four- velocity ti*, and 
its reference counterpart ^^^^ are availed of. Despite its four-dimensional 
clothing, equation (|1.8| ) defines an effectively three-dimensional stress tensor: 
since both the auxiliary metric and its reference countepart are orthogonal 
to the four-velocity, one finds 

(1.10) T.ijU^ = 0, 

as it is appropriate for describing the elastic stress. Carter and Quintana 
extended Rayner's theory to a non-linear regime [0], as it is necessary when 
dealing with astrophysical situations. Nordstrom derived |p his equations 
of motion from Hamilton's principle and showed their coincidence with the 
ones obtained through the general relativistic law of conservation for the 
energy tensor 

(1.11) T;f = 0. 

The later authors mentioned above and many others [^]-|17] instead relied 



only on equation ( 1.11| ) for setting up their equations of motion. There is 



however a notable exception in the Lagrangian formulation given |18| by 
Kijowski and Magli. 

2. A REALLY FOUR-DIMENSIONAL EXTENSION OF HOOKE'S LAW 

The diverse proposals for a general relativistic reformulation of Hooke's 
law all adopt an effectively three-dimensional stress-strain relation. It seems 
therefore natural to wonder whether a truly four-dimensional Hooke's law 
can have some physical sense. The general relativistic extension of equa- 



tion (1.6) is formally immediate: in a pseudo-Riemannian space-time whose 
metric is gik one considers a contravariant four- vector field ^'^{x ) that aims 
at representing some "displacement", and uses its covariant counterpart to 
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define a four-dimensional "deformation" tensor 

(2.1) Sik = -{ii;k + ik;i)- 

A four-dimensional "stiffness" tensor density C*'^'™ is then introduced; it 
will be symmetric in both the first pair and the second pair of indices, since 
it will be used for producing a "stress-momentum-energy" tensor density 

(2.2) T^'^ = a^^'^Sim. 

through the four-dimensional extension of Hooke's law. If T is meant 
to be the overall energy tensor density, it must obey the four-dimensional 



conservation law ( 1.11 ) , that of course has to do with the equations of motion 
of the field ^'^{x^). From a formal standpoint, this is nearly the end of 
the story. The real question is: can this abstract scheme find a physical 
interpretation, at least in some limit condition? 

Let us start inquiring whether the vector field ^^{x^) can admit of the 
following physical meaning. One considers a co-ordinate system such that, 
at a given event: 

(2.3) gik='nik = diag{l,l,l,-l), 

while the Christoffel symbols are all vanishing, and the components of the 
four-velocity of matter are 

(2.4) u^ =u'^ = u^ = 0, u^ = 1. 

In this co-ordinate system we suppose it possible to measure, at the cho- 
sen event, both the components of the spatial displacement of the elastic 
medium from some relaxed condition and the proper time, that we interpret 
as temporal displacement. If these four quantities can be determined at any 
event with this procedure, they can be used to define, through the appropri- 
ate transformations, the four- vector field ^^{x ) in an arbitrary co-ordinate 
system. Let us assume that the field ^*(a;'^) that we have previously intro- 
duced in a purely formal way really allow for this physical interpretation. 
Then, if the material is unstrained in the ordinary sense, the only nonva- 
nishing component of Sik should be 544, whose value should be —1. This 
remark suggests defining the four- velocity n* of matter through the equation 

(2.5) e,kn' = u\ 

A necessary condition for the above definition to hold is: 
(2.6) detiC^k - SI) = 0. 

This shall be one equation that the field ^* must satisfy; the number of 
independent components of ^* will thereby be reduced to three. We observe 
that the very definition of a four- velocity field n'(x ) by starting from the 
field ^* is not ensured a priori; its possible existence will depend on the 
physical properties of the model that one considers. 



A four-dimensional "stiffness" tensor C^^''"^ possibly endowed with phys- 
ical meaning can be built as follows. We assume that in the locally Min- 
kowskian rest frame defined above the only nonvanishing components of 
(jikim gj.^ (jXuar -^[^^^ ^]-^g tentative role of elastic moduli, and 



(2.7) C*^^^ = -p, 

where p measures the rest density of matter (energy) . We need to write the 
four-dimensional "stiffness" tensor in an arbitrary co-ordinate system. The 
task can be easily accomplished if matter is isotropic when looked at in the 
above mentioned Minkowskian rest frame. Just for the sake of simplicity, we 
shall deal henceforth only with this case. One avails of the auxiliary metric 
( |1.9D ; then the part of (7*'^''" stemming from the ordinary elasticity of the 



isotropic medium can read [15] 



(2.8) Cf/™ = -AfS'"" - /"(/7^™ + ^"'t*^'), 

where A and p are assumed to be constants. The part of (7*^'™ accounting 
for the inertial term reads instead 

(2.9) C^'"^ = -pu'u''u^u"'. 

When no other fields are present, the tensor density T*'^ defined by equation 
r2[ ) must fulfil the four conditions ( 1.11| ). Since u* is defined in terms of (,^ 



and of the metric tensor go. via equation (p.5|), imposing the four conditions 
( |1.11| ) means postulating four nonlinear partial differential equations that, 
for a given gn-, must be obeyed by the three independent components of the 
field ^* and by the scalar field p. Thanks to ( p.8[ ) and (p.5|), the putative 
elastic part of the energy tensor T*^ for isotropic matter reads: 

T^i. = c'^'^^Sim = -Kg'"' + u'u'^Wn^ - 1) 

(2.10) -Pie"' + t'' + uiiu't^'' + n'^e''')], 

and is orthogonal to the four- velocity, as an intrinsically three-dimensional 
tensor should be, while the supposedly inertial part of T*'^ turns out to be 
effectively so. In fact, due to equation (|2.5|): 

(2.11) Tt = C'^^'Si^ = pu'uK 



When A = /x = we have no elasticity at all, and the field equations ( |1.11 ) 
reduce to: 

(2.12) Wu^};fc = {pu^};fcu* + pu'.^u^ = 0. 
By contracting the last equation with Uj one finds: 

(2.13) {pu^],k = 0. 

This equation is convenient for finding p, once the problem of motion, en- 
tirely written in kinematic terms: 

(2.14) u\f,u^ = 0, 
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has been solved. The world lines of matter are then timelike geodesies of 
the spacetime whose metric is gi^. Once the geodesic equations are solved 
for u* also the displacement field ^* is in principle determined through the 



equations ( |2.5| ) and (2.6). In the limit case when equation ( |2.3| ) holds ev- 
erywhere the world lines of matter are straight and by performing a Lorentz 
transformation we can put matter everywhere to rest in the co-ordinate sys- 
tem, say, X* = {x,y,z,t). Then equation ( |2.4D is satisfied everywhere; we 
can fulfil both (|2.5D and (p. 61) by posing: 



(2.15) e = iO,0,0,t). 

3. The classical theory of elasticity as limit case 

Let us hold fixed both the co-ordinate system and the metric selected 
above, and assume that now the constants A and ;U are nonvanishing, but 
that M* differs very slightly from its rest form ( |2.4D . We deal with the 
components u^ as with first-order infinitesimal quantities, while of course 
the increment to u^ shall be infinitesimal at second order. Since one expects 
that ^p, like ^'^, and ^^ will be at most first order infinitesimal quantities, 
the first order approximation to equation (|2.5| ) requires 

(3.1) u'' = i% 
and 

(3.2) u^ = i% = 1. 

From the above remark about ^*p and from equations (|3.l[) , (|3.2|) one gathers 
that ^\ will fulfil equation (|2.6D to the required first order. Let us check 
the changes occurred to equation ( 2.13| ), that expressed the conservation of 



matter when ordinary elasticity was absent. One writes: 



(3.3) 






T%u^ -- 


= (r'=u,);fc-r^n 


i;k 


and, since 










(3.4) 








T^jtuk = 0, 




equation 


(3.3|) 


comes 


to read 






(3.5) 






(- 


pu%k-Tlui.k = 


0. 



While the inertial term of this equation is a first-order infinitesimal, the 
putative elastic term contains only infinitesimal quantities of a higher order. 
Therefore equation ( 2.13| ) still holds, with the required approximation. We 



write now the field equations ( 1.11 ) by retaining only the first-order terms. 
The equations for i = 1,2,3 are: 

{pu'^u^ - Xirj'^'' + u''u''){^';;^ - 1) 

(3.6) -^i[C'' + C"'" + uK^^'e'-' + t^'e''")]},fc = 0. 



Due to ( |2.13| ) and ( |2.5D , one can rewrite these equations as 

.pi 



(3.7) pu-^^u^ = XC^^ + ^{e,'^P + e,''^-) 



one eventually obtains, to the required order: 

(3.8) <4,4 = Ae"x + Mr" + e'''n.p- 



Equation (3^) says that ^4 can differ from unity at most for second order in 

:P> 

':P 



finitesimal terms; hence, with the required accuracy, ^'^f can be substituted 



for ^,m in the previous equations, that come to read: 

(3.9) P^^^ = K'':'; + ^^{e^'' + in,p■ 

The field equation for i = 4 reads: 

(3.10) -^,[e'^ + e^^ + uKA''" + ^'e''')]},fe = 0. 

The elastic terms provide only second-order contributions; therefore only 
the inertial term survives and, to the required order, the last field equation 
reads 

(3.11) {pu\^],k = 0, 



as it occurs for the geodesic motion. The conservation equation ( 2.13| ) could 



have been adopted as fourth equation as well. The three equations 



and the conservation law ( 2.13 ) exactly match the corresponding equations 



of the classical theory of elasticity. 

4. Conclusion 

As far as one can understand from the results of the previous Section, 
postulating a truly four-dimensional Hooke's law does not seems to be a 
false step, since it allows to merge both linear elasticity and inertia into a 
sort of extended, four-dimensional elasticity. Investigating the content of 
this approach when the metric is highly nonflat and the relative speed of 
different portions of matter is comparable with the velocity of light is a 
mathematically difficult and physically useless undertaking, since one must 
expect that the linear approximation of the classical theory of elasticity will 
fail to account for the behaviour of matter under such extreme conditions. 
Writing the equations of motion of elastic matter in presence of a weak, wavy 
deviation oi gik from fiatness|jis instead both mathematically affordable and 
slightly more useful: in the complete lack of experimental evidence about 
the way elastic matter interacts with gravitational waves, it will at least 
allow for a comparison with the theoretical predictions achievable through 
the diverse approaches recalled in the Introduction. 



^The very concept of weak gravitational wave is not without pitfalls. As wittily stressed 
by Eddington |19| , one can easily manifacture perturbations of the flat metric that "prop- 
agate" with the speed of thought. 
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Confronted with the bareness of the experimental landscape, one can 
however seek consolation in the sense of intellectual fulfilment that looking 
at known things from a new perspective always provides. The expectation 
aroused by the analogy between the four-dimensional constitutive equation 
of electromagnetism ( |1.5[ ) and the three-dimensional expression ( |1.6D of the 
stress-strain relation of classical elasticity was not void. Viewing inertia as 
a kind of elasticity is indeed possible, and this possibility is not constrained 
to the weak field, slow motion approximation. The reformulation of the 
inertial term of T*^ was possible thanks to equations ( |2.5D and ( |2.6D , whose 
role was central in forcing the four- vector field ,^* to account for both inertia 
and ordinary, classical elasticity. 
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